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We study beating of a probe field with a time-varying susceptibility in a coherently prepared 
Raman medium. We consider the general case of an arbitrary variation of susceptibility, which 
corresponds to a superposition of an arbitrary number of excited Raman transitions. We derive a 
general analytical solution and conservation relations for this process. We show that the interference 
between Raman polarizations may substantially affect frequency modulation and pulse compression 
for the probe field. 
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We have recently presented a detailed analysis of ul- 
trashort laser pulse compression by parametric beating 
with a sinusoidal molecular oscillation jl). In this Brief 
Report, we extend the previous work to the general case 
of an arbitrarily complex molecular motion. We assume 
a wave-like molecular excitation, produced by external 
fields in a dispersionless medium. Molecular motion re- 
sults in a time-varying susceptibility, which (in the ab- 
sence of dispersion) leads to time variations in phase and 
group velocities for a probe field Q. We analyze the ef- 
fect of these variations on pulse propagation, derive an 
analytical solution in terms of an integral for a time trans- 
formation, and use it to obtain conservation relations. 

It was shown recently that coherently driven molecular 
oscillations can produce frequency modulation with ul- 
trabroad bandwidth (2[ , and result in subfemtosecond 
jl], ||, ||] and subcycle j7| pulse compression. A molec- 
ular oscillation can be either sinusoidal (a single coher- 
ently excited Raman transition) ^ H, ^, (| or more 
complex (multimode vibrational or rotational molecular 
wave packets, which correspond to a coherent phased ex- 
citation of many Raman transitions) ||, [|. Molecular 
wave packets can be excited either impulsively |], ||, |), 
or by applying several quasi-monochromatic laser fields 
0. Generalization of our previous results (l) to the im- 
portant case of a complex molecular motion will provide 
insight into the Raman pulse compression techniques. 

The essence of our result is simple: A pulse is com- 
pressed or stretched by time-varying phase and group 
velocities, such that envelope compression or stretching 
is accompanied by frequency increase or decrease, re- 
spectively. Pulse amplitude and oscillation frequency are 



changed by reciprocal factors, such that the photon num- 
ber, the pulse area, and the number of optical oscillations 
are conserved, while the pulse energy is not. When the re- 
sultant change in the pulse bandwidth is small compared 
to the optical carrier frequency, pulse deformation by 
the complex non-sinusoidal molecular motion can be ob- 
tained as a multiplication of the effects of individual har- 
monic components of the motion. However, for stronger 
modulation this intuitive thinking can not be applied. In 
order to emphasize this point, in the second part of our 
Report we consider an example of a biharmonic molecu- 
lar excitation and show that frequency modulation and 
pulse compression are significantly affected by the inter- 
ference between two Raman polarizations. 

We study the propagation of a probe field E in a Ra- 
man medium, characterized by a polarization P. In the 
local time r = t — z/c, the reduced wave equation reads 
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We assume that dispersion is negligible. The instanta- 
neous susceptibility is defined as x = P/cqE. We as- 
sume that X = Xo + Xmij — z /v), where xo is the linear 
time-independent susceptibility, Xm is the Raman time- 
varying susceptibility, and v — 2c/ xo is the average phase 
and group velocity in the local time coordinates. 

We use the reduced local time r\ = r — z/v. In the 
coordinates z and r], Eq. ([!]) becomes 

I^^-^XMEM. (2) 

The solution to Eq. (|^) is found to be 

E(z,r]) = E ia (s)G(ri), (3) 
where E{ n (s) = E(z = 0, s) is the input field, s is the 
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input time determined from the output time rj by 
rr > d9 z 



(4) 



Xm{6) 2c' 
and G(?7) is the compression factor given by 

Xrniv) 

Differentiating Eq. (Q) with respect to rj, we obtain 

According to Eqs. (||) and (0), the height and duration of 
optical oscillations in the probe field are changed by the 
reciprocal factors E(z,rj) / ' E- m (s) — G{rf] and dr}/ds = 
1/G(r]), respectively. The value G(rj) > 1 (G(rj) < 1) 
indicates pulse compression (stretching) in the vicinity 
of rj. From Eqs. (||) and @), we find the relation 
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E(z,r))dr) = / E in (s)ds, 

>)i J Si 
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which describes the conservation of the pulse area. 

We find from Eq. (|J) that the instantaneous oscillation 
frequency of the output field is 



w osc (?y) = G(r])ui , 



(8) 



where luq is the input frequency As seen, u; osc (77) is mod- 
ulated in time by the compression factor G(ri). Using 
Eqs. d), (|), and (|), we find 



2 



Tiu> osc (ri) 
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This relation describes the conservation of the photon 
number, which is always satisfied for Raman processes. 

We introduce the mean frequency uj{j]\,T}2) = (j]2 ~ 
?/i ) _1 ui sc(v) drj for the oscillations in the time inter- 
val (771 , 772) • Then, we find from Eqs. (||) and (|J) that 
w(?7i, 772) (?72 — rji) = u>q(s2 — si), that is, the product of 
the pulse length and the mean frequency is constant dur- 
ing the propagation process. Furthermore, Eq. (|J) says 
that E(z,ri) = if Ei D (s) — and vice versa, that is, a 
zero of the input field at an input time s leads to a zero 
of the output field at the corresponding output time 77. 
Hence, the number of optical oscillations is conserved. 

To get insight into the behavior of the compression 
factor G, we derive an explicit, approximate expression 
for this factor. For this purpose, we consider a particular 
case where the deviation of an output time 77^ from its 
corresponding input time Si is small compared to the 
characteristic variation time T m of Xmi that is, r)i\ <C 
T m . We take rj close to rji, and use the approximation 
Xm(9) = Xm(Vi) + x' m (Vi){0-Vi) to calculate the integral 
in Eq. (g). Then, we find 
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Hence, Eq. d§) yields 
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Equation (If]) shows that the factor G is approximately 
determined by the time derivative of the Raman suscep- 
tibility Xmi multiplied with the propagation length z. 
Note that, around the chosen time rji, the factor G is 
multiplicative, that is, pulse deformation by the complex 
molecular motion can be obtained as a multiplication of 
the effects of individual components of the motion. The 
assumption |sj — rji\ <C T m requires 



exp [-{z/2c) X ' m {m)] ~ 1 
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Condition ( [f2] ) is satisfied when (z / c)\x' m {r]i)\ "C 1 and 
(z I ' c)\xm{Vi)\ ^ ^m, i.e., when the Raman susceptibil- 
ity and its modulation are small. In this case, we have 
G(r)i) = 1 — (z/2c)x' m ( r li) — 1) indicating that the pulse 
compression is weak. Condition ( |f2| ) is also satisfied 
when r\i is a zero of Xm ■ 

In general, G is not multiplicative with respect to 
the individual susceptibilities. To show this, we con- 
sider the case where the Raman susceptibility Xm is 
a sum of two sinusoids corresponding to two individ- 
ual Raman transitions, namely, x-m = Xa + Xb, where 
Xj = xf^ sin^-Ty + ipj) for j — a,b. The compression 
factor corresponding to the individual component Xj is 
given by the periodic comb function 0, [Ic[] 



G, 



e «^cos 2 ^^ +e- 
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where otj = (ujj/2c)Xj ■ In general, the total compres- 
sion factor G is different from the product of the indi- 
vidual compression factors G a and Gb- The reason is 
that the compression factor G is produced by the time- 
varying Raman susceptibility Xm via a nonlinear mech- 
anism, see Eqs. (||) and (|J). Due to the nonlinearity of 
this mechanism, the individual susceptibility components 
may interfere with each other in beating with the field. 

To see the interference between the individual suscepti- 
bility components in beating with the field, we illustrate 
the functions G and G a Gb in Figs. g(a) and 0(b), re- 
spectively. The plots are calculated for the parameters 
a a z — 0.8, ctbz = 0.6, w& = 0.07 u> a , and ip a = ipb = 0. 
The time is normalized to the Raman period T a — 2-k /ui a . 
Both G and G a Gb reveal a sequence of teeth correspond- 
ing to the oscillation of the susceptibility component Xa, 
which has a higher frequency. The peaks of these teeth 
are modulated in accordance with the oscillation of the 
component Xb, which has a lower frequency. The differ- 
ences between G and G a Gb are clearly observed. 

To see in detail the difference between G and G a Gb, 
we magnify in Fig. || the curves of Fig. |f| for the time 
windows f < ?7/T a < 2 (a) and 7.5 < r]/T a < 8.5 (b). We 
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FIG. 1: Comparison between the total compression factor G 

(a) and the product of individual compression factors G a Gb 

(b) . The plots are calculated for the parameters a a z — 0.8, 
othZ = 0.6, uib = 0.07 uj a , and ip a = <fb = 0. 
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FIG. 2: Magnification of the curves for G and G a Gb of Fig. |l| 
in the regions 1 < r)/T a < 2 (a) and 7.5 <rj/T a < 8.5 (b). The 
temporal profiles of a short probe field at the output in these 
time regions are also shown. The carrier frequency and pulse 
length of the input field are u>o = 15.2 cu a and T = 0.08 T a , 
respectively. The input peak time is rj p = 0.83 T a (a) or 8.8 T a 
(b). The peak value of the input field is one unit. 



also plot in Fig. |2| the output field for the cases where 
the input peak time is rj p = 0.83 T a (a) and 8.8 T a (b). 
The carrier frequency and pulse length of the input probe 
field are ujq = 15.2 ui a and T — 0.08 T a , respectively. We 
observe G < 1 < G a G b and G > 1 > G a G b around the 
times r)/T a — 1.5 and 8.1, respectively. These relations 
imply that, due to the interference between \a and Xb, 
the total compression factor G and the product of in- 
dividual factors G a Gh may give opposite indications on 
pulse stretching and compression. Thus, the interference 
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FIG. 3: Comparison between the mixture (a), the cascade 
M a Mb (b), and the cascade MbM a (c) of the individual Ra- 
man systems M a and Mb- The parameters for the media are 
a a z = otaZa = 0.8, a.bz = (XbZb = 0.6, ujb = 0.07 uj a , and 
fa = <fib = 0. The carrier frequency and pulse length of the 
input field are the same as for Fig. H. The input peak time is 
rj p = 3.1 T a - The peak value of the input field is one unit. 



between the susceptibility components can substantially 
affect frequency modulation and pulse compression. 

In the above, the total compression factor G has been 
compared with the product G a G b of the individual com- 
pression factors. However, G a G b cannot be rigorously 
considered as a compression factor for the case where 
there is no interference between Xa and Xb- Below we 
will compare the compression factor of a combined sys- 
tem with the compression factor of a sequence of two 
individual systems. 

We call M a Mb a cascade in which the probe field is 
sent through, first, a cell M a (with susceptibility Xa and 
length z a ) and, then, a cell M b (with susceptibility Xb 
and length z b ). For this configuration, we have E ou t(rf) = 
E\ a {s)G ab (rj) and ds/dr] = G ab {r]), where 

G ab (r 1 )=G a (r 1 ',z a )G b (r 1 ,z b ). (14) 

Here rj and s are determined by the equations |l], |l0) 

tan[(u;6?/ + <p b )/2] = e- abZb tan[(w fc ?7 + ip b )/2], 
tan[(w a s + ip a )/2] = e- aaZa tan[(cj a ?/ + ip a )/2], (15) 

Note that G ab is not symmetric with respect to the in- 
dices a and 6, i.e., G ab ^ G ba . In other words, the cas- 
cades M a M b and M b M a produce different compression 
factors. This feature is different from the results of Rcf. 
0, which are valid for the case of limited bandwidths. 

Comparison between Eqs. (||) and ( |l4| ) shows that the 
compression factor G of the combined system is different 
from the factors G ab and G ba of the cascades. To see 
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the difference between the mixture and cascades of the 
individual systems M a and Mf,, we plot the factors G, 
G ab , and G ba in Figs. ||(a), |(b), |(c), respectively, for 
the time window 4 < rj/T a < 5. We also plot the output 
field and the output spectrum for the input peak time 
rj p = 3.1 T a . The other parameters are the same as for 
Fig. g. We observe G a b < G = 1 < Gba around the 
time rj/T a = 4.5. These relations as well as the temporal 
and spectral profiles of the output field show clearly the 
differences between the combined and cascade systems in 
frequency modulation and pulse compression. 

In the case where ctjZj <C 1 for j = a, 6, we find to 
lowest order in cxjZj the compression factor G a b = 1 — 
J2j a j z j cos(ujjri+ipj) and the input-output time relation 
s = j] — *}2j{otjZj /uij) sm(u) jT] + Hence, with the use 

of the generating function sxn 6 — YlkL-oa Jk(Qe lke of 
the Bessel functions we can expand an input oscilla- 
tion e luJ ° s into a series of output harmonics exp[i(o;o + 
qjU)j)r)]. Then, we can approximate the spectrum of 
the field at the output of a series of M a and Mb as the 
product of the Bessel-function spectra of the individual 
cells 0. The condition otjZj <C 1 means that the pulse 
bandwidths jjZjUj |J, with jj = ajLUo/uij, produced 
by the individual cells, are small compared to the optical 
carrier frequency luq. The same results are also obtained 
for the case of mixed cells (except that individual cell 
lengths Zj should be replaced by a common length z). 
Under the condition of small pulse bandwidths, there is 
no interference between \a and %{,, and therefore, no dif- 
ference between the mixture and cascades of M a and Mb 
in frequency modulation and pulse compression. 

Finally, we demonstrate a numerical example for a real- 
istic system, namely, for a cell containing mixed H2 and 
D2 molecules. We take ou a = 587 cm" 1 and ojb = 179 
cm" 1 so as to correspond to the rotational transitions 
of H2 and D2 . We use a probe with a carrier frequency 
wo = 20 000 cm -1 , a pulse length T = 4.5 ps, and a peak 
time rj p = 0. We assume the comb depths a a z — 0.587 
and abz = 0.179, which correspond to the modulation 
depths 7 a z = jbZ = 20. Similar to the results of Ref. 
for the case of a cascade, a broad spectrum with a large 
number of sidebands, ranging from 8000 cm" 1 to 46 000 
cm" 1 , is generated. Unlike the results of Ref. |7j], the 
spectrum in Fig. ^(a) is asymmetric. Such asymmetry 
is due to the deviation of the spectrum from the Bessel- 
function spectrum, and occurs when the comb depths 
ctjZ are not too small compared to unity. The output 
field in Fig. |(b) is a long train of slightly compressed 
sections. As has been shown in Ref. |7j for the case of a 
cascade, single-cycle pulses can be synthesized when the 



sidebands are phase-corrected by a phase compensator. 
We illustrate in Fig. |](c) a single-cycle pulse synthesized 
from the phase-corrected spectrum of Fig. ||(a) . 

In summary, we have studied beating of a probe field 
with a time-varying susceptibility of a Raman medium. 
We have derived a general analytical solution and conser- 
vation relations for this process. We have shown that the 
interference between Raman polarizations may substan- 
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FIG. 4: Spectrum (a), output field (b), and time-domain syn- 
thesis by a phase compensator (c) for a cell containing H2 
and D2. The parameters for the medium are uj a = 587 cm -1 , 
LOb — 179 cm" 1 , a a z = 0.587, atz — 0.179, and <p a = ft = 0. 
The carrier frequency, pulse length, and peak time of the in- 
put field are ljq — 20 000 cm" 1 , T — 4.5 ps, and r\ v — 0, re- 
spectively. The peak value of the input field is one unit. The 
synthesized field is obtained when the sidebands are phase- 
corrected to the same phase of zero. 



tially affect frequency modulation and pulse compression. 
We emphasize that the analysis of this work can be di- 
rectly applied only when dispersion is negligible. When 
dispersion becomes substantial, our analytical solution 
and conservation relations still provide a useful insight, 
but a numerical simulation of pulse propagation is nec- 
essary in order to obtain an exact solution. 

We thank A. K. Patnaik for helpful discussions. A. V. 
S. acknowledges support from the Welch Foundation. 
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